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Abstract — We apply an approach for estimating Value -at-Risk (VaR) describing the tail of the conditional distribution of a 
heteroscedastic financial return series. The method combines quasi-maximum-likelihood fitting of AR(1)-GARCH(1,1) model 
to estimate the current mean as well as volatility , and extreme value theory (EVT) to estimate the tail of the adjusted 
standardized return series. We employ the approach to investigate the existence and significance of the calendar anomalies: 
seasonal effect and day-of- the-week effect in Americas Indexes VaR. We also examine the statistical properties and made a 
comprehensive set of diagnostic checks on the one decade of considered Americas Indexes returns. Our results suggest that 
the lowest VaR of considered Americas Indexes negative log returns occurs on the fourth season among all seasons. 
Moreover, comparatively low Wednesday VaR is captured among all weekdays during the test period. 

Keywords — Risk Measures , Value -at-Risk, GARCH models, Extreme value theory, Generalized Pareto Distribution, Day- 
of-the-week effect, Seasonal effect. 


I. Introduction 

In today’s financial world, the large increase in the number of traded assets in the portfolio of most financial institutions has 
made the measurement of market risk a primary concern for regulators and for internal risk control. Following the Basle 
Accord on Market Risk (1996) most banks in more than 100 countries around the world have to calculate its risk exposure 
for every individual trading desk. Banks are also required to hold a certain amount of capital as a cushion against adverse 
market movements. Value-at-Risk has become the benchmark risk measure. From mathematics point of view, VaR is simply 
a quantile of the profit-and-Loss distribution of a given portfolio over a prescribed holding period. The importance of VaR is 
undoubted since regulators accept this model as a basis for setting capital requirements for market risk exposure. 

In this paper, we discover the calendar anomalies in Americas equity market movements, which including the seasonal effect 
and the day-of-the-week effect on Americas Indexes returns. The calendar effect in stock market returns includes day-of-the- 
week effect, weekend effect, January effect, and holiday effect, etc. It has been widely studied and investigated in finance 
literature. Studies by Cross (1973), and Rogalski (1984) demonstrate that there are differences in distribution of stock returns 
for each day of the week. Studies by Baillie and DeGennaro (1990), Berument and Kiymaz (2001) posit that day-of-the-week 
effect has an impact on stock market volatility. In recent years, another stream of research has considered seasonality in stock 
returns and volatility, see Saunders (1993), Bouman and Jacobsen (2002), Hirshleifer and Shumway (2003), Kamstra, 
Kramer and Levi (2003), and Cao and Wei (2005), etc. These studies generally report that calendar anomalies are present in 
both returns and volatility equations in the stock market. None of these studies, however, test for the possible existence of 
day-of-the-week and seasonal variation in stock return VaR. Hence, the goal of this paper is to characterize the VaR of 
Americas Indexes returns. Based on investigations of the day-of-the-week effect and seasonal effect in extreme risk, we also 
provide valuable and applicable analysis for equity market investors. The major obstacle to this investigation is a viable 
measure of tail risk over time. 

We are concerned with tail estimation for those considered financial return series. Our basic assumption, whose validation is 
examined in this paper, is that returns follow a stationary time series model with stochastic volatility structure. The presence 
of stochastic volatility implies that returns might dependent over time. Therefore, we consider to model the return 
distribution as the conditional return distribution where the conditioning is on the current volatility and mean. Although VaR 
only characterizes the extreme quantiles, disregarding the centre of the distribution, estimation of the extreme quantile is not 
an easy task. As one wants to make inferences about the extremal behavior of a portfolio, there is only a very small amount 
of data in the tail area of a sample set. Advanced methods and tools are needed to enable us to explore beyond the range of 
the limited data set. In this study, we use a semi -parametric method, based on extreme value theory (EVT), which is rather 
effective for obtaining reliable estimates. 

EVT is widely used to study the distribution of extreme realizations of a given distribution function, or stochastic processes 
that satisfy suitable assumptions. The foundations of the theory were laid by Fisher and Tippett (1928) and Gnedenko (1943), 
who demonstrated that the distributions of the extreme values of an independent and identically distributed sample from a 
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cumulative distribution function, when adequately rescaled, can converge towards one out of only three possible 
distributions. Unfortunately, most financial time series are not independent, but exhibit some very delicate temporal 
dependence structure. In this study, we capture it by a fully parametric method, which is based on an econometric model for 
volatility dynamics and the assumption of conditional normality, AR-GARCH model. We use AR(1)-GARCH(1,1) model 
and quasi-maximum-likelihood estimation to obtain estimates of the conditional mean and the conditional volatility. 
Statistical tests and exploratory data analysis confirm that the standardized returns, i.e. mean and volatility adjusted returns, 
do form approximately i.i.d. series. If we only use GARCH model to estimate VaR , the assumption of conditional normality 
does not seem to hold for real data. Thereafter, we use threshold methods from EVT to estimate the distribution of the 
standardized returns. EVT is a well known technique in many fields of applied sciences including risk management, 
insurance and engineering. Numerous research studies surfaced recently which analyze the extremes in the financial markets 
due to currency crises, stock market turmoils and credit defaults. The behavior of financial series tail distributions has, 
among others, been discussed in Reiss and Thomas (1997), McNeil and Frey (2000), Longin (1999) and (2000), and Ameli 
and Malekifar (2014). An estimate of the conditional return distribution is now easily constructed from the estimates of the 
conditional mean and volatility as well as the estimated distribution of the standardized returns. We learned the central idea 
of the dynamic two stage extreme value process from McNeil and Frey (2000), to forecast daily VaR with historical data in a 
moving window. This approach reflects two stylized facts exhibited by most financial return series, stochastic volatility, and 
the non-normal behavior of conditional return distributions. 

In terms of the organization of this paper, we first introduce the data set and present a comprehensive set of diagnostic checks 
on it in Section 2. We review certain aspects of VaR and introduce the extreme value based approach to it in Section 3. 
Section 4 contains the empirical tests and results on the Americas Indexes returns. The seasonal effect and day -of-the- week 
effect on Americas Indexes returns VaR and economic implications of the empirical results are examined in Section 4 as 
well. Conclusion is made in Section 5. 

II. Data Exploration and Statistical Analysis 

2.1 Data Description 

Our sample covers the period from July 17, 2006 to November 13, 2015. Five different Americas Indexes, namely, the S&P 
500, Financial Select Sector SPDR ETF, NASDAQ-100 Technology Sector, Dow Jones Utility Average, and Dow Jones 
Transportation Average, are used to characterize the performance of specific sectors of the market. The S&P 500 is an 
American national index composed of large capitalization stocks. It represents the overall performance of the stock market. 
The Financial Select Sector SPDR ETF tracks the overall S&P Financial Select Sector Index. The NASDAQ-100 
Technology Sector is an equally weighted index based on the securities of the NASDAQ -100 Index that are classified as 
Technology according to the Industry Classification Benchmark classification system. The Dow Jones Utility Average is a 
stock index from Dow Jones Indexes that keeps track of the performance of 15 prominent utility companies. The Dow Jones 
Transportation Average is a U.S. stock market index from S&P Dow Jones Indices of the transportation sector, and is the 
most widely recognized gauge of the American transportation sector. The collections of those indices’ daily adjusted closing 
price were from Yahoo Finance. The adjusted closing price is used to develop an accurate track record of the stock’s 
performance. 

Furthermore, we use the daily negative log return to examine extreme losses of the stock. Let P t denotes the adjusted closing 
price of a stock on day t , then the daily percentage change on the day is defined by 

r, = -lOOtojffj, = lOOto* — . 

”■ (i) 

The reason for using the negative returns is that we are mainly interested in the possibility of large losses rather than large 
gains. 
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Figure 1: Time plots of the Standards and Poors index from 2006-7-17 to 2015-11-13. There are 
2351 DAYS EXCERPT FROM THE SERIES OF ADJUSTED CLOSING PRICE (UPPER) AND THE SERIES OF NEGATIVE 

DAILY LOG RETURNS (LOWER). 

Fig.l shows the time plots of adjusted closing price and negative daily log returns of S&P 500 from July 17, 2006 to 
November 13, 2015. From the lower plot, we observe that daily log returns of the Index show clear evidence of volatility 
clustering. That is, periods of large returns are clustered and distinct from periods of small returns, which are also clustered. 
If we measure such volatility in terms of variance, then it is nature to think that variance changes with time, reflecting the 
clusters of large and small returns. We also observe that there are more pronounced peaks than one would expect from 
Gaussian data. Since the possibility of time-varying variance and non-normal behavior are noticed in Fig.l, we provide 
formal tests to check the stationarity, normality, and independency of those log return series. 

2.2 Test for Stationary Property 

The invariance of statistical properties of the return time series corresponds to the stationarity hypothesis that the joint 
probability distribution of the returns does not change when shifted in time. Here we use the KPSS test to verify the 
hypothesis of weak stationarity, i.e. time invariance of the mean value and the autocorrelation function of America Indexes 
return series. 

Proceeding in the spirit of Kwiatkowski, Phillips, Schmidt and Shin (1992), we assume that the series 1 can be 
decomposed into the sum of a random walk and a stationary error 1. We express this symbolically by writing 

r, = CT a +- 


1 In general, the assumption of KPSS test is that the series {rt} can be decomposed into the sum of a deterministic trend, a 
random walk and a stationary error. In this study, we only consider to test series stationary with no trend. 
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where is assumed to be stationary; at is a random walk, i.e. ^ 1 " Here ^ 1 Ms a white noise series with zero 

mean and variance ^ . The hypothesis for the KPSS test is 




i.e. the series is stationary vs 




.a,; * 0 


, i.e. the series is not stationary 


j 

S. 1,2, 

= r — Gt F 

Assume 1 1 1 as residuals of the regression of ' on an intercept, as the partial sum 

process of the residuals. The KPSS test statistics is 


KPSS = 



( 2 ) 


] |^y 

where ^ is a consistent estimator of the long-run variance of 1 . 

The rejection rule is that if the value of the KPSS statistic in Eq.(2) exceeds the critical values estimated in [19], or the p- 

rt H 

value is less than or equal to the significance level , we reject 11 . 


Table 1 

Stationary Test and Normality test results on the five Americas Indexes negative daily log 

RETURNS FROM 2006-7-17 TO 2015-11-13 


Data 

S&P 

XLF 

NDXT 

DJU 

DJT 

Observations 

(2350) 

(2350) 

(2398) 

(2351) 

(2350) 

KPSS Test for time 

series level stationarity 




KPSS statistics 

0.17128 

0.28586 

0.072364 

0.081087 

0.1073 

p - value 

0.1 

0.1 

0.1 

0.1 

0.1 

Shapiro- Wilk Test for time series 

normality 




W statistics 

0.87934 

0.82416 

0.94139 

0.88601 

0.9479 

p - Value 

<2.2e-16 

<2.2e-16 

<2.2e-16 

<2.2e-16 

<2.2e-16 


H ct 

n . In this study, we choose the significant level ^ as 5% for KPSS test. 

The KPSS rest results on the five America Indexes negative daily log returns from July 17, 2006 to November 13, 2015 
are shown in Table 1, all p-values are greater than the significant level 5%. Therefore, we accept the null hypothesis and 
conclude that the five return series are stationary during the test period. 

2.3 Tests for Normality 

In studying the financial time series, people usually assume that the log return follow a normal distribution. However, it is 
unrealistic to make the normality assumption on America Indexes returns, according to our QQ-plot on the S&P 500 negative 
daily log returns against the normal distribution, see Fig. 2. 
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Figure 2: Quantile-Quantile plot of the S&P 500 negative daily log returns from 2006-7-17 to 

2015-11-13 AGAINST THE NORMAL DISTRIBUTION 

To prove that the returns do not follow normal distribution, we use one of the most powerful formal normality tests (as 
Razali and Wah (2011) demonstrate) - the Shapiro- Wilk test, to verify an empirical fact that the five America Indexes 
negative daily log return series do not have the normality property. 


The Shapiro-Wilk test utilizes the null hypothesis principle to check whether the series 
distributed population. The Shapiro-Wilk test statistic is defined as 




come from a normally 


<2 



2,^-7.? 


r-1 


jr tt 

where r is the i-th order statistic; f is the sample mean; 


(A, 


are the weights. 


(3) 


The value of W lies between zero and one. Small values of W lead to the rejection of normality whereas a value of 1 
indicates the normality of data. Consequently, we reject the null hypothesis if the p -value of the test is less than the 
predetermined significance level, which is 5% in this study. 

Applying the Shapiro-Wilk test on the five America Indexes negative daily log returns from July 17, 2006 to November 13, 
2015, we show the test result in Table 1. Because all p-values are less than 2.2e-16, we reject the null hypothesis and 
conclude that all the five return series are not normally distributed during the test time period. 


2.4 Test for Correlations 


Besides the verified stylized fact of the fat tail distribution, we further explore the autocorrelations correlations for the returns 

U l 

series as well as their squared values. We begin by considering the autocorrelation function of a time series 1 s 1 . The 
correlation between "■ and its past values ^ 1 is called the lag-1 autocorrelation of ^ and is commonly denoted by ^ ,! , 
for 1=0,1,. . .t. Under the weakly stationary assumption, we assume P ■' is a function of 1 only. 


P: = 


Cov(r r ) _ Co v(f r , r t _. ) y, 


,.) Van.*) y* 


(4) 
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where the property } 7fi f or a weakly stationary series is used. 

R 


Jp l T 

For a given sample of returns 1 1 : 1 , let 
can be represented as: 


}T 

be the sample mean. The lag-1 sample autocorrelation of 1 1 1 s 1 


P t = 


(a) 


ilK-r?UT- 1 } 


-J}<t<T- I. 



(b) 


( 5 ) 



Figure 3: Sample autocorrelations of (a) returns and (b) squared returns of the S&P 500 

(2006-7-17 TO 2015-11-13) 

If a time series is not autocorrelated, the estimate of ^ will not be significantly different from 0. 

Fig. 3(a) shows that the sample autocorrelation coefficient ^ plotted against different lags 1 (measured in days), along with 
the classical 95% significance bands around zero for S&P 500 negative daily log returns. The dashed lines represent the 

, 1 .96 

upper and lower 95% confidence bands , where the time length for our S&P 500 returns is ^ ” 2350 days. A 

stylized fact that absence of autocorrelation for the daily price variations is illustrated in Fig. 3(a). The series of S&P 500 
returns displays small autocorrelations, making it close to a white noise. However, the S&P 500 squared returns are strongly 
autocorrelated, see Fig. 3(b). This property is not incompatible with the white noise assumption for America Indexes returns, 
but shows that the white noise is not strong. 


Based on the statistical analysis for the five America Indexes negative daily log return series, we discovered that those 
America Indexes returns are stationary, and uncorrelated time series, yet are not normally distributed and the squared returns 
are strongly correlated. Those properties illustrate the difficulty of daily price returns modeling. Any satisfactory statistical 
model for daily returns must be able to capture the main stylized facts, including the leptokurticity, the unpredictability of 
returns, the existence of positive autocorrelations in the squared returns, and the conditional heteroscedasticity. Some 

Jj - V 

computations of VaR are based on the assumption that the series 1 : 1 is normally distributed, or has t-distribution, see 
reference [27] [1] [3] [18] [15]. That is the main reason why these study can use volatility to estimate VaR. However, the real 

time series may not follow any known distributions. To overcome the difficulty of a return series with an 

unknown distribution, we compute the VaR of America Indexes returns by the Extreme Value Theory, which avoid making 

assumption about the distribution of ^ ^ . 
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III. Value-at-Risk Methodology 

Exposure to risk can be summarized as a single number by estimating the VaR, which is defined by Jorion as “the worst 
expected loss over a great horizon within a given confidence level”, it is crucial to have an accurate estimate on VaR. Before 
we go into our VaR methodology in more detail, we introduce the probability theory behind VaR. 

3.1 VaR Introduction 

VaR is the amount that might be lost in a portfolio of assets over a specified time period with a specified small failure 
probability a , usually set as 0.01 or 0.05. In this paper, we choose such short period as one day. Suppose a random variable 
X characterizes the distribution of negative daily returns of a portfolio, the right-tail a-quantile of the portfolio is then defined 

to be the ^ : !j ^' ' such that 


Fr{X<VaRJ=\-a. 


( 6 ) 


VaR 1 _ rt 

The ^ is the largest value for X such that the probability of a loss over a day is no more than L u - . Consequently, the 

can be viewed as a point estimate of potential financial loss. 

Although the parameter ^ is arbitrarily chosen, the analysis in this study does not refer to the process of choosing the 
parameter of VaR which were considered to be C {0.0! , 0.05 * 0.1 j Hence the crux of being able to provide an 
accurate estimate for the VaR j s [ n estimating the cutoff return ^ : !j ^' ' . 

Following the approach by Longin (1999a,b), and McNeil and Frey (2000), we use a two-stage approach to estimate the VaR 
of the five America Indexes negative daily log return time series. 


(1) Fit a AR(1)-GARCH(1,1) model to the returns and use a quasi-maximum-likelihood approach to estimate parameters. 
Use the fitted model to standardize the raw returns to a strict white noise process, i.e. independent, identically distributed 
process with zero mean and unit variance. 


(2) Use EVT to model the tail of the marginal distribution of the standardized returns, and use this EVT model to estimate 

VaR . 


3.2 Standardization - Estimating and ^ 1 by QML 

Since stock returns have heavy-tailed and/or outlier-prone probability distributions, we use GARCH models to deal with both 
the conditional heteroskedasticity and the heavy-tailed distributions of American Indexes returns, we consider the reason for 
outliers may be that the conditional variance is not constant, and the outliers occur when the variance is large. According to 
our test, the returns process is not Gaussian. Therefore, we assume the standardized, i.e. mean and variance adjusted 
American Index returns series is an i.i.d. white noise process with a generalized Pareto distribution. 


Let L 1 1 be a strictly stationary time series representing the negative daily log return on a financial asset price. We fix 

a constant memory n so that at the end of day T our data consist of the last n negative daily log returns ^ L ' ' ' ' ^ 1 . 


Let ^ be the information about the return process available up to time t. We assume that the dynamics of 
a realization from a AR(1)-GARCH(1,1) process, which are given by 


{r,}L 


to be 


r. = u r ■+ a r z, > 
ft = ^ -l- $r ( , 
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where the innovations JLi are white noise process with zero mean, unit variance, and marginal distribution function F; 

W> 0„Ot > 0, an( j P > (J . ^ conditional mean ^ 1 , and the conditional volatility ^ ^ 1 are measurable, ^ 1 1 
the information about the return process available up to time ^ — 1 . 


is 


This model is fitted using the quasi-maximum-likelihood estimation (QML) method, which assumes normal distribution 
and uses robust standard errors for inference. It means that the likelihood for a GARCH(1,1) model with normal innovations 

is maximized to obtain parameter estimates ^ ^ ^ ^ . Although this amounts to fitting a model using a distributional 
assumption we do not necessarily believe, the QML method delivers reasonable parameter estimates. Bollerslev and 
Wooldridge (1992) proved that if the mean and the volatility equations are correctly specified, the QML estimates are 
consistent and asymptotically normally distributed. From Eq.(7), we get estimates of the conditional mean 

r ' 1 ■ ‘ " ^’ r ^ and the conditional 


volatility 


A .-A A} ofUlL 


. The estimated conditional volatility of the S&P 500 


returns derived from the GARCH fit 2 is shown in Fig. 4. 



Figure 4: Estimate of the conditional standard deviation derived from QML fitting of AR(1)- 
GARCH(1,1) MODEL OF THE S&P 500 NEGATIVE DAILY LOG RETURNS (2006-7-17 TO 2015-11-13). 

To check the adequacy of the model and to use in next stage of the approach, we calculate the standardized returns 



°V- N+l (8) 


(a) 


(b) 


- 

r ■ r; :: : 

c 

.2 ~ 
s 0 

a; 

fc r 
O n 
O 

..4. 1 1 . ... ' ... 

3 n 


< 


S 


<M 




Lag 


Lag 


Figure 5: Sample autocorrelations of (a) standardized returns and (b) squared standardized 

RETURNS OF THE S&P 500 (2006-7-17 TO 2015-11-13). 


2 

In this research, the AR(1)-GARCH(1,1) model is fit to the American Indexes negative daily log returns using R’s garchFit 
function in the fGarch package[32]. 
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The standardized returns should be i.i.d. if the fitted model is tenable. In Fig. 5, we plot the sample autocorrelation of the 
standardized S&P 500 returns as well as the squared standardized S&P 500 returns. As shown in Fig. 3, while the raw returns 
are clearly not i.i.d., this assumption may be tenable for the standardized returns. 

We end the standardization stage by calculating estimates of the conditional mean and variance for day ^ + 1 , which are the 
1-step forecasts 


A.-i = ft +■ i 


3.3 VaR Estimation - Apply Extreme Value Approach 

In the second stage, we estimate the upper tail behavior of the cumulative distribution function of the standardized returns F 
by using the Extreme Value Theory (EVT). EVT is experiencing a boom in the financial field, especially with respect to the 
application to the market risk measure VaR. Its appearance as a popular instrument for estimating VaR can be explained as a 
consequence of two factors. On the one hand, the assumption of normality of financial markets does not reflect the reality of 
the situation. As a consequence, the VaR estimation methods which based on the normality assumption may provide 
inaccurate estimates. Historical or Monte Carlo simulation methods arise as alternative methods. But given the difficulties 
and “slowness” of these methods, EVT has been used as a powerful solution for estimation of VaR. On the other hand, 
although VaR can be calculated with simulation methods, it still has limitations, so this measure needs to be complemented 
with others. We present it as a way of solving the problem of fat tails when calculating VaR. 


The mathematical foundation of EVT is the class of extreme value limit theories, originally stated by Fisher and Tippett 
(1928) and later derived rigorously by Gnedenko (1943). The central result in EVT is that the extreme tail of a wide range of 
distributions can approximately be described by the Generalized Pareto distribution (GPD), which is derived by Smith(1989), 
Davison and Smith (1990). 


For a random variable X, we first fix some high threshold ! and consider the distribution of excess values 


Y=X-p 


as 


F£y) = PiU- p<y\x>ji)= ffcZl ^ 

i -Fiji) 


( 10 ) 


where F is the underlying distribution of X, ^ is the conditional excess distribution function. Pickands (1975) introduced 
the GPD as a two parameter family of distributions for exceedances over a threshold. More precisely, it was proved that for a 


large class of underlying distribution functions F, the conditional excess distribution function 




as 


fi — > oiy = : F(x) cl} 


, is well approximated by 




H a 

where 


Ay) 


is called GPD, specified as 


N 


/V}= I -{! + £ — ) 

tj 



( 11 ) 


The parameters of GPD are the scale parameter and the shape parameter 


EVT describes specifically at the distribution of the standardized returns in the tails. The tail fatness of the distribution is 


reflected by the shape parameter: the case when 


6 


<0 


means thin tails, 


£=0 


means the kurtosis is 3 as for the standard 


normal distribution; while ^ ^ ^ implies fat tails. Therefore, the shape parameter measures the speed with which the 
distribution’s tail approaches zero. The fatter the tail, the slower the speed and the higher the shape parameter is. Since 
almost all returns in EVT assume that the returns are i.i.d., the analysis was developed on the standardized returns which, in 
many cases, could be reasonably assumed to be i.i.d.. Because we are interested in extreme negative returns, we use EVT to 
model the right tail of the distribution, i.e. the standardized returns in excess of a high threshold. 
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It is necessary to choose a specific threshold to confine the estimation to those observations that are above the given 
threshold. However, it is difficult to apply threshold based methods because of the lack of a clear-cut criterion for choosing 
the threshold. If the threshold is chosen too low, the GPD may not be a good fit to the excesses over the threshold, and 
consequently there will be a bias in the estimates. Conversely, if the threshold is too high, then there are not enough 
exceedances over the threshold to obtain reliable estimates of the extreme value parameters, and consequently, the variances 
of the estimators will be high. In this paper, an optimal threshold is selected by employing graphical methods, known as the 

Hill plot and the mean excess plot. The Hill plot displays the estimated values of shape parameter ^ as a function of the cut- 
off threshold in order to find some interval of candidate cut-off points that yields stable estimates of the shape parameter ^ . 
Technical details about Hill plot can be found in Hill (1975). The mean excess function is the mean of exceedances over a 
threshold. If the underlying distribution of those exceedances follows a GPD, then the corresponding mean excess must be 
linear in the threshold. Details about the mean excess plot are described in Davison and Smith (1990). Fig. 6 shows the Hill 
plot and the mean excess plot of the negative daily S&P 500 log returns. A threshold 1.968748, with 132 exceedances, seems 
to be reasonable for the S&P 500 returns. 


Threshold 

4d30 1.973 1.270 0.903 0.B32 



Order Statistics 


Mean Excess Plot 



Figure 6: Hill plot (left) and Mean Excess Plot (right) of the S&P 500 negative daily log 

RETURNS (2006-7-17 TO 2015-11-13) 


r c? 

We have seen that the GPD contains two parameters, shape parameter ' and location parameter u . They can be estimated 
by using either parametric or non-parametric methods. From the research of Hosking and Wallis (1987), for the tail index 

^ ^ , it can be shown that maximum likelihood regularity conditions are fulfilled and that maximum likelihood estimates 


^ } based on a sample of n excesses are asymptotically normally distributed. Therefore, we use the parametric 

approach, maximum likelihood method (MLE), to estimate parameters in GPD. 


Next, we make explicit the relationship between excess value and the standardized return series, denoted as ^ . We may 

use the following relationship to estimate the VaR of the standardized asset returns ^ . 

Assume that ^ ,! are i.i.d. random variable with CDF F, and a high enough threshold ^ is given. Define 

jV. = card{t i}. Then 




_ } F{ti+y)~F(u) 

1 -Fiji) 


F ;i {y) = l'r(z i ~fi>y\z ( 

i.e. 


>fi)= 


i >o 

Fip) 


which is equivalently to 
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FUi + y) = F{u)F(y). 


Then, the estimators of ^ and ^ ^ can be written as: 


n “ r? 

>{>■)= I- tf,,, (>■) = {1 + | ) t "% 


( 12 ) 


r £7 r C? 

where ^ and " are maximum likelihood estimators of the shape parameter and the location parameter " . Therefore 

the tail estimator can be written as: 


a (7., 


(13) 




This relationship between probabilities allows us to obtain VaR for the original asset return series 1 ; 1 . More precisely, for a 
specified small probability K such that the — ^ upper tail quantile VaR of ^0 j s V . Consequently, for a given small 
probability ® , one can check that the VaR of holding a long position in the asset underlying return I 


is 


FT * ') 

a=lM>> |>v)=FrU r ,^ V )= t^ni ~ P* ' — — — > fi) = F(— — -JJ) = — Sli — ? 

°V,] °7«] 17 tii F W 


VaR, = { 


a., 

V-+-? 
L C 


(T 








, n m ) 


(14) 


We favor the extreme value approach, or the GPD approach in this study to tail estimation mainly for three reasons. One is 
that in finite samples of the order of points from typical return distributions, EVT quantile estimators are more efficient than 
the historical simulation method. Moreover, considering the fact that most financial returns series are asymmetric, the EVT 
approach is advantageous over models which assume symmetric distributions such as t-distributions, GARCH distribution 

family. In addition, comparing with Hill method which is designed specifically for the heavy tail ( * ^ ^ ) data, the EVT 

approach to VaR has larger applicability since it also applicable to light tail ( ^ ^ ) cases or even short tail ( ^ ^ ^ ) cases. 

IV. Empirical Results and Calendar Effect Analysis 


We backtest the approach on the five Americas Indexes historical series of negative daily log returns: the the Standard and 
Poors index S&P 500, the financial index SPDR ETF, the technology index NASDAQ- 100, the utility index Dow Jones 
Utility Average, and the transportation index Dow Jones Transportation Average. As introduced in Section 2, we excerpt all 
five indexes’ adjusted closing price from July 17, 2006 to November 13, 2015. 


To backtest the two stage approach, we first estimate and ^ and use it to standardize the daily negative log returns 

^ ^ . The reason we use the negative returns is that loss occurs when the returns are negative for a long financial position. 
We show the estimation results of AR(1)-GARCH(1,1) model for Americas Indexes negative daily log returns in Table 2. 


After getting the standardized returns 
returns. 




, we apply the second stage to test the calendar effect on Americas Indexes 
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Table 2 

Estimation Results of AR(1)-GARCH(1,1) model for Americas Indexes negative daily log 

RETURNS FROM 2006-7-17 TO 2015-11-13 


Parameter 

S&P 

XLF 

NDXT 

DJU 

DJT 

Mean equation AR(1) 
00 

-0.069493 

-0.037231 

-0.082928 

-0.042878 

-0.060223 

0 

-0.058278 

-0.079456 

-0.015441 

-0.029874 

-0.016125 

Variance equation G ARCH (1.1) 





vj 

0.025183 

0.027637 

0.026652 

0.018436 

0.028343 

a 

0.114753 

0.128632 

0.078668 

0.097303 

0.077237 

d 

0.867333 

0.867966 

0.908884 

0.889307 

0.910656 

LogLikelihood 

-3342.299 

-4189.893 

-4062.424 

-3268.324 

-4117.66 

jh Value of Standardised 

Residuals Tests 




Shapiro-Wilk Tfcst, 

0 

1.92972e-15 

6.5971 62e-ll 

9.209965e-13 

7.049808e-lC 

Ljung-Box Test Q(10) 

0.1749702 

0.02211829 

0.4162687 

0.6369799 

0.3831675 

Ljimg-Box Test Q(15) 

0.1254892 

0.03006419 

0.4009677 

0.8425757 

0.2452127 

Lj ung- Box Test Q(20) 

0.2380792 

0.06046306 

0.4091626 

0.9131188 

0.1832375 

Ljimg-Box Test Q(10) 
for squared residuals 

0.02556951 

0.4572453 

0.7371797 

0.06534186 

0.08447315 

Ljung-Box Test Q(15) 
for squared residuals 

0.05896544 

0.4860042 

0.8815682 

0.1705194 

0.2434914 

Ljung-Box Test Q(20) 
for squared residuals 

0.1596955 

0.5549062 

0.7223037 

0.3516701 

0.2387464 


4.1 Seasonal Effect on Americas Indexes VaR 


Because a three-month period on a financial calendar acts as a basis for the reporting of stock earnings and the paying of 
dividends, the seasonal effect is a vital factor in determining stock performance. To identify the existence of seasonal effect 

on Americas Indexes returns, we divide each Americas Index’s standardized returns ' into the following four subsets: 

{z.licQK r= L2J.4, 


which are referred to the four quarters Index returns. A quarter refers to one-fourth of a year and is typically expressed as 
Table 3 provides a summary of descriptive statistics for those considered return series. 

Table 3 

Descriptive Summary Statistics of Americas Indexes Seasonal Standardized Returns 


Summary Statistics 

S&P 

XLF 

NDXT 

DJU 

DJT 

Observations n 

Q1 

550 

550 

565 

550 

550 

Q2 

569 

569 

579 

569 

569 

Q3 

620 

626 

639 

626 

626 

Q4 

604 

604 

614 

605 

604 

Skewness 

Qi 

0.81 

0.6 

0.45 

0.42 

0,32 

Q2 

0.56 

0.34 

0.16 

0.45 

0.4 

Q3 

0,53 

0,27 

0.35 

0.35 

0,21 

Q4 

023 

0,03 

0.16 

0.19 

0,24 

Kurtosis 

Q1 

3,15 

2,97 

1.19 

1,3 

1.23 

Q2 

0,56 

0,68 

0.33 

0*54 

0,47 

Q3 

1,32 

1*5 

0.95 

1.7 

0,59 

Q4 

0.73 

0.87 

0.39 

0.35 

0.49 

Shapiro-Wilk normality test p- Value 
Ql 1.021e-10 1.959e-09 

2.88e-06 

3.952e-06 

3.669e-05 

Q2 

1.37e-07 

0.0004483 

0.001367 

1.454e-05 

0.0001729 

Q3 

3.42e-09 

L923&07 

3.879e-05 

2.109e-06 

0.003301 

Q4 

0.0001126 

0.005276 

0.06336 

0.006695 

0.02605 


Page | 163 


International Journal of Engineering Research & Science (IJOER) 


[Vol-1, Issue-9, December- 2015] 


Table 3 reports skewness and kurtosis for the standardized return series of each quarter. In statistics, skewness and kurtosis, 
which are normalized third and fourth central moments of a process, are often used to summarize the extent of asymmetry 
and tail thickness. For the normal distribution, kurtosis is 3. We observe that distributions of all four seasons’ standardized 

returns are positively skewed, indicating that they are nonsymmetric. Further, except for the S&P 500 ^ standardized 
returns, all kurtosis are less than 3, indicating that those series have distributions with tails that are thinner than those of the 
normal distribution. This indication of non-normality is also supported by the Shapiro-Wilk test results, which reject the null 
hypothesis of a normal distribution at 5% significance level. Based on the Ljung-Box test results from Table 2 and the 
Shapiro-Wilk test results from Table 3, we consider all five Americas Indexes seasonal standardized returns are i.i.d. and 
non-normally distributed processes. 


Table 4 

Results from Fitted GPD for Standardized Returns & Estimates for VaRs for Negative 

Daily Log Returns 


Standardized Return 
(Threshold /i) 

R&P 

{1.651907) 

XL.F 

(1.658118) 

MDXT 

(1.595358) 

D.IU 

(1.630336) 

D.TT 

(1,590404) 

Slltf.pfl PiirfiTTH’tnT £ 

Q1 

0.01263226 

0.08324439 0.2507820 

-0.2261632 

-0.05614722 

Q2 

-0.4155444 

0.04520885 -0.8514770 

-0.2794200 

-0.06843322 

Q3 

-0441,8723 

-0.1682080 -0.0387090 

0.1329915 

-0.2223172 

Q4 

-0 1534931 

0.1375718 

-0.1974187 

0.00310129 

-0.1744704 


Scale Parameter a fi 


Qi 

0,77437170 

0.77008703 

0.4372119 

0.9268009 

0,68909844 

Qa 

0,9013507 

0.535 71 577 

0.9237556 

0.7847136 

0.01659084 

Q3 

1.2207484 

0.8877799 

0.6161069 

0.6354024 

0.7442377 

Q4 

0.6340782 

0.4532491 

0.6423728 

0.49605948 

0.6551974 

Exceedances j 

Qi 

37 

36 

44 

31 

39 

Q2 

44 

46 

45 

44 

43 

Q3 

41 

44 

48 

35 

46 

Q4 

34 

33 

35 

33 

32 

05% quantile VaR (VaRu.uo) of original negative daily log returns 
QI 1 1.61934 1,947586 2.00777 1.996732 

1,008087 

Q2 

1.739842 

2.000061 

2.166308 

2.252225 

1.921836 

Q3 

1.694324 

2.036937 

2.059414 

1.951028 

1.947205 

Q4 

1.474158 

tc 

IQ 
f- 
1 — ■ 

1.864537 

1.917086 

1.691029 

99% quantile VaR (YaJQ. (l , ) of original negative 
QI 1 2.797445 3.453345 

daily log returns 
3.147227 3.450343 

3,032853 

Q2 

2.562514 

3,002238 

2.820178 

3.334415 

2.009932 

Q3 

2.857024 

3,333057 

3.172462 

3.341999 

2,047104 

Q4 

2.301884 

2.665381 

2.879703 

2.932556 

2.675514 


Thereafter, we apply the extreme value approach to the considered seasonal standardized return series. Table 4 summarizes 

r £7 

estimation results of the shape parameter ^ , the scale parameter l ' from fitted GPD model for the standardized Americas 
Indexes returns as well as 0.95 quantile VaR and 0.99 quantile VaR for the original considered negative daily log returns. To 
better investigate systematic seasonal differences for Americas Indexes VaR, we first find a proper threshold for each 
seasonal series by Hill plot and then choose the highest one as the common threshold for all four seasonal series. Using the 
same threshold allows a better comparison of the quarterly VaRs. 

From Table 4, we find that for S&P 500, at a quantile level of 95%, the smallest estimated VaR among all seasons is 

04 CM 

1.474158 for the “ returns; at a quantile level of 99%, the smallest estimated VaR is 2.301884 for the ^ returns as well. 
This implies that, in the fourth season, with the AR(1)-GARCH(1,1)-GPD model, we are 95% confident that the expected 
overall Americas equity market value would not lose more than 1.474158% for the worst case scenario; we are 99% 
confident that the expected market value of the S&P 500 would not lose more than 2.301884%. Similar interpretations can be 
made for the other Americas Indexes. 
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In comparison of all four seasonal returns, it is also interesting to note that our model produced the smallest VaR in the fourth 
season, at the 95% quantile level for all five America Indexes. While at the 99% quantile level, except for the NASDAQ -100 
technology index, the four seasonal VaRs exhibits analogous characteristics as observed from 95% quantile VaRs under 
different Americas Indexes seasonal returns. Moreover, given the quantile levels, the corresponding VaR estimates for S&P 
500 seasonal returns are less than the rest Indexes seasonal returns. It indicates that the trading risk of S&P 500 is the 
smallest among all five Americas Indexes. 

Our findings have important implications for investors and financial institutions. For example, for conservative investors 
who would prefer lower risk, they can choose to trade during the lower VaR period or trade lower risk stocks to avoid 
potential high loss. 

4.2 Day-of- the- Week Effect on Americas Indexes VaR 

To formally test the timing and existence of weekly patterns, we divide the whole standardized returns J to five subsets 
by day-of- the- week, which is written as: 


{z : \d t =]}J = WAJ, 


where are dummy variables such that if day t is a Monday £t,s ^ , remove the data if ‘^ M ^ ; if day t is a Tuesday 

^ , otherwise remove the data, etc. The five subsets are the considered Indexes weekly returns for Monday through 
Friday respectively. The basic statistical characteristics of the five return series are calculated and shown in Table 5. 


Table 5 

Descriptive Summery Statistics of Americas Indexes Weekly Standardized Returnes 


Summary Statistics 

S&P 

XLF 

XDXT 

DJU 

DJT 

Observations n 

Monday 

442 

442 

468 

442 

442 

Tuesday 

480 

480 

483 

481 

480 

Wednesday 

483 

483 

485 

433 

483 

Thursday 

473 

473 

479 

473 

473 

Friday 

471 

471 

482 

471 

471 

Skewness 

Monday 

0.41 

0.24 

0,22 

0.17 

0.41 

Tuesday 

0.77 

0-68 

0.31 

0.47 

0.12 

Wednesday 

0.04 

0.27 

0.08 

-0.1 

0.08 

Thursday 

0.18 

0.03 

0.38 

0.47 

0.47 

Friday 

0-65 

0.43 

0.49 

0.83 

0.53 

Kurtosis 

Monday 

1.22 

1.15 

1.28 

2.04 

0.8 

Tuesday 

3.1 

3.14 

0,88 

0.62 

0.3 

Wednesday 

0.47 

1.01 

0.61 

0.41 

0.57 

Thursday 

0.94 

0.78 

0.44 

0.77 

0.52 

Friday 

1.04 

1.1 

0,41 

1.82 

1.69 

Shapiro- Wilk normality test p - Value 
Monday 2.07c-06 0.0008976 

0,0004619 

1. 911o-05 

0.0007413 

Tuesday 

2.242e-09 

1 .94M19 

0.001296 

9.745e-05 

0.5958 

Wednesday 

0.002344 

0.001305 

0,005957 

0,4103 

0.00485 

Thursday 

3. 49 5c- 08 

0.006136 

0,0004381 

8. 0660 -O 5 

0.0001246 

Friday 

8J6Le-08 

8.494e-05 

6.554e-06 

2.965e-09 

4.802e-(J6 


Except the Dow Jones Utility Average Wednesday standardized returns, the distribution of the rest weekday’s standardized 
returns are slightly positively skewed, indicating that they are nonsymmetric. The kurtosis of the S&P 500 and the SPDR 
ETF Tuesday standardized returns show Gaussian property while the rest of the week returns are substantially departure from 
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normal distribution. The Shapiro-Wilk test results also indicates that normal distribution is not a realistic assumption for the 
weekly standardized returns for considered Indexes. 

Next, we apply the VaR estimation approach to capture the day-of-the-week effect on Americas Indexes VaR and show the 
result in Table 6. Given Americas Indexes and the 99% quantile level, we captured the comparatively low risk in 
Wednesday. Among the five Americas Indexes, the number of exceedances is comparatively small in Monday and Friday. 


V. Conclusion 
Table 6 

Results from Fitted GPD for Standardized Returns & Estimates for VaRs of Negative Daily 

Log Returns 


Standardized Return 
(Threshold 

S&P 

(1.685503) 

XT.F 

(1.597324) 

KDXT 

(1,616329) 

d,tt; 

(1,498996) 

DJT 

(1.654929) 

Shape Parameter £ 
Monday 

-0.490609 

0.08324459 

0.5182840 

0.2425261 

-0,4065323 

Tuesday 

0.09756212 

0.1224898 

0.0999698 

-0,1381534 

-0.1377256 

Wednesday 

-0.1117621 

0.1075106 

-0.2249917 

-0.1726002 

-0.2138769 

Thursday 

-0-3322608 

-0.1846136 

-0.3968165 

0.06310120 

-0.1077049 

Friday 

-0,5091703 

-0,1816384 

-0,1592645 

-0.174639 

-0.09100887 


Scale Parameter \ 7 . ± 


Monday 

1.239719 

0.77968703 

0.5634457 

0.4779636 

1.0962016 

Tuesday 

0.65606741 

0.7608726 

0.5040591 

0.6698328 

0.6560630 

Wednesday 

0.5055943 

0.4589546 

0.6706803 

0.5467823 

0.6060613 

Thursday 

-0.3322698 

0.7758208 

0.9214631 

0.49605048 

0.6201783 

Friday 

1. 1695172 

0.8252605 

0.5825002 

1,037172 

0.87551725 


Exceedances A 1 ,., 


Monday 

j- 

25 

30 

20 

26 

22 

Tuesday 

53 

56 

53 

37 

27 

Wednesday 

50 

55 

27 

44 

40 

Thursday 

38 

30 

30 

33 

35 

Friday 

24 

28 

30 

32 

17 

95% quantile 

VaR (VaUo.u&) uf original negative daily log returns 


Monday 

1.1574378 

1.9127 

1.8863157 

1.844478 

1.714618 

Tuesday 

1 .63393 

1.094884 

1.976933 

2.085929 

1.802911 

Wednesday 

1.537115 

1.838576 

1.875154 

2,123752 

2,034037 

Thursday 

1.827274 

2.054521 

2.272032 

1,949974 

1,981222 

Friday 

1 .4.56647 

1.802248 

2.10879 

2.121393 

1 .406727 

99%: quantile 

VaR (Valuta) uf original negative daily lug returns 


Monday 

2.784977 

3.425015 

2.848478 

5.020582 

3,117782 

Tuesday 

2.733362 

3.5,54609 

3-044567 

3.180491 

2.810845 

Wednesday 

2.214732 

2.7,58804 

2.914875 

2.9,56329 

2.836568 

Thursday 

2.776619 

3.134498 

3,535698 

2,981338 

2,945308 

Friday 

2.642806 

3.025975 

5.009405 

3.77907 

2,867233 


With the empirical analysis of this paper we demonstrate how we can use a GARCH based dynamic EVT approach to model 
VaR for short term forecasting. The dynamic EVT method has the advantage of dynamically reacting to changing market 
conditions which is useful in getting better VaR forecasts. We apply the two stage approach on five Americas Indexes 
negative daily log return series. Empirical findings in this paper show that both seasonal effect and day-of-the-week effect 

04 

are present in Americas Indexes returns. We captured the comparatively low VaR in “ and Wednesday for considered 
returns during the test period. 
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Overall, our findings have implications for investors, financial institutions, and futures exchanges. For example, for 

conservative investors who would prefer lower risk, they can choose to trade during the lower VaR period to avoid potential 

high loss. The two stage approach to VaR can also be used in other stock or asset returns. Finally, it has significant value for 

investors and regulators in terms of an in depth analysis of the equity market. 
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